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Abstract. We present a full treatment of the microcanonical ensemble of the ideal hadron-resonance gas 
in a quantum-mechanical framework which is appropriate for the statistical model of hadronization. By 
using a suitable transition operator for hadronization we are able to recover the results of the statistical 
theory, particularly the expressions of the rates of different channels. Explicit formulae are obtained for the 
phase space volume or density of states of the ideal relativistic gas in quantum statistics which, for large 
volumes, turn to a cluster decomposition whose terms beyond the leading one account for Bose-Einstein 
and Fermi-Dirac correlations. The problem of the computation of the microcanonical ensemble and its 
comparison with the canonical one, which will be the main subject of a forthcoming paper, is addressed. 

PACS. 12.40.Ee - 05.30.-d 



1 Introduction 

The revived interest in the statistical model of hadron pro- 
duction is mainly owing to its application to heavy ion col- 
lision where an equilibrated source of hadrons is expected. 
This model has given strikingly good results in elementary 
collisions as well and this finding has triggered some 
debate about their interpretation [2]. A proposed one is 
that hadronization occurs at some critical energy density 
[3] (or maybe another related parameter) of a number of 
massive pre-hadronic colourless extended objects (hence- 
forth referred to as clusters) which are formed as a result 
of the underlying non-perturbative strong-interaction dy- 
namics and which thereafter decay coherently into mul- 
tihadronic states 0]. In this scheme, the single cluster's 
decay rate into any channel would be determined only by 
its phase space with no special dynamical weight (phase 
space dominance). Thereby, the observed statistical equi- 
librium would not be the effect of a collisional thermaliza- 
tion process between formed hadrons over long-lived ex- 
tended regions in the final state, rather of equal quantum 
transition probabilities from a cluster to all accessible final 
states. By accessible it is meant that one must comprise 
only those states fulfilling conservation laws, i.e. having 
the same quantum numbers as the initial cluster's. The set 
of states with fixed energy-momentum, angular momenta, 
parity and internal charges is denned as microcanonical 
ensemble, though the same name is usually employed to 
denote the set of states with fixed energy-momentum and 
internal charges, relaxing angular momentum and parity 
conservation. We will not make any distinction either; 
whether the constraints of angular momenta and parity 
are meant to be included, it will be clear from the con- 
text. 

Although the microcanonical ensemble is the correct 
statistical ensemble to use in hadronizing a single cluster, 



so far all actual data analyses within the statistical model 
have been carried out in the framework of the canoni- 
cal or grand-canonical ensemble of the hadron gas, i.e. 
with hadronizing sources described in terms of a temper- 
ature and taking into account the conservation of energy- 
momentum and angular momentum only on average. The 
microcanonical ensemble has been used very seldom , 
mainly owing to the hard and long computations involved. 
Indeed, in high energy collisions, where many clusters are 
produced, the use of the canonical ensemble is favoured 
by fluctuations of masses and volumes, which tend to re- 
duce the importance of exact conservation of energy and 
momentum. It is even possible that fluctuations make the 
system of many clusters equivalent (as far as Lorentz in- 
variant quantities, such as average multiplicities, are con- 
cerned) to a large global cluster obtained by ideally clump- 
ing them |2j . While the canonical ensemble is in fact a bet- 
ter and better approximation of the microcanonical one 
for large values of cluster's mass and volume, we have no 
quantitative estimate of how large they ought to be 1 , so 
the use of the canonical ensemble is, until now, justified by 
the agreement with the data. On the other hand, it would 
be desirable to have a more precise and quantitative as- 
sessment of the goodness of the canonical approximation. 
An explicit calculation of the microcanonical ensemble of 
the hadron gas is also necessary if we want to test the 
statistical model at lower energies (say y/s < 10 GeV) 
where conservation laws are expected to play a major role 
and canonical approximation is not a good one. Further- 
more, it would be very useful having at our disposal a 
Monte-Carlo algorithm for microcanonical hadronization 
of single clusters in high energy collisions to be used for nu- 
merical calculations of quantities for which an analytical 



1 Recently, a calculation has been done for pp collisions with 
a restricted set of hadrons .6: 
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expression cannot be obtained. Thus, providing a reliable 
and fast numerical algorithm for the calculation of micro- 
canonical ensemble of the hadron gas and comparing the 
results with the canonical approximations are the main 
goals of this work which will be described in two papers. 

In this first paper, we will confine ourselves to the 
analytical development of the microcanonical formalism, 
while the numerical calculations will be the main subject 
of the second paper. In fact, another major motivation of 
this work is to provide a consistent formulation of the sta- 
tistical hadronization model starting from quantum tran- 
sition probabilities, which is still lacking despite the ap- 
parent simplicity of the picture and the fact that its foun- 
dations were laid more than 50 years ago [Hj ■ The attempts 
to derive the statistical theory results from S'-matrix under- 
suitable hypotheses, were mainly carried out by Hagedorn 
[§] and Cerulus rTU| on the basis of time-reversal argu- 
ments. However, the reasoning is quite involved in this 
approach and one needs separate treatment of dynamical 
matrix element averaging for multiplicities and spectra. 
The inclusion of quantum statistical effects in the micro- 
canonical ensemble of the relativistic hadron gas has been 
done consistently only for large volumes and not suffi- 
ciently general. In fact, in this traditional treatment, par- 
ticle states confined in the cluster's volume are assumed to 
be eigenstates of energy-momentum, which is true only if 
the volume is so large that the entailed energy-momentum 
uncertainty can be neglected, what is not generally the 
case when dealing with small volumes (this will be dis- 
cussed more in detail at the end of Sect. 2). Furthermore, 
in that approach, the whole treatment did not start from 
the statistical theory of multiple production and it is thus 
not easy to generalize if angular momentum and parity 
conservation are to be included. 

Therefore, we believe that a coherent general reformu- 
lation of the statistical model of hadronization is needed. 
In this paper, we will go along the whole formalism start- 
ing from the basic assumptions and will recover some well 
known formulae in literature, like those in ref. fTI], as ap- 
proximations of more general ones in case of sufficiently 
large volumes. In particular, we will recover the TV-body 
relativistic phase space expression (without angular mo- 
mentum and parity conservation) without treating con- 
fined particle states as energy-momentum eigenstates, an 
assumption which is correct only asymptotically. We will 
show how conservation laws are to be implemented in the 
most general case, thus providing an usable framework 
to obtain more general expression of TV-body phase space 
when angular momentum and other conserved quantities 
are to be taken into account. Furthermore, we will ex- 
plicitely show how the microcanonical ensemble reduces 
to the canonical one for large cluster's volume and mass. 
The formulae presented in this work will be then the basis 
of the numerical computations in the second paper |12|. 

The paper is organized as follows: in Sect. 2 we will 
introduce a basic formulation of the statistical hadroniza- 
tion model and follow the path leading to the microcanon- 
ical ensemble; in Sect. 3 we will develop in detail the mi- 
crocanonical formalism for an ideal hadron-resonance gas 



with full quantum statistics; in Sect. 4 the microcanonical 
partition function will be calculated and the approxima- 
tions needed to obtain closed expressions stressed, while 
in Sect. 5 the transition from the microcanonical to the 
canonical ensemble described; finally, in Sect. 6 the calcu- 
lation of the physical observables will be discussed. 



2 Statistical hadronization of a cluster 

The fundamental assumption of the statistical hadroniza- 
tion model is that the final stage of a high energy colli- 
sion results in the formation of a set of extended colour- 
less massive objects, the clusters or fireballs, producing 
hadrons in a purely statistical manner: that is, all multi- 
hadronic states within the cluster volume and compatible 
with cluster's quantum numbers are equally likely. Clus- 
ters can indeed be thought as very short-lived extended 
resonances, much alike to bags of the bag model |13| . 
They differ from clusters proposed in other hadronization 
models [H] as they are endowed with a spacial exten- 
sion. In this picture, the cluster's decay rate into a given 
iV-particle channel should be proportional to the number 
of multiparticle states within the volume V of the clus- 
ter, which can be expressed, in the limit of Boltzmann 
statistics and neglecting angular momentum and parity 
conservation, as: 

N 

rocV N /d 3 Pl ...dV<5 4 0P-E^) C 1 ) 

J i=i 

where P is the four-momentum of the cluster. The dis- 
tinctive feature of the statistical model is essentially the 
appearance of a finite volume in the decay rate, which 
makes the above expression different from that of the de- 
cay rate of a massive particle usually found in textbooks, 
where asymptotic states are defined over an infinitely large 
volume: 

i— 1 

While these two equations have in principle a different 
physical meaning, there have been several attempts to de- 
rive an equation like (JIJ from Eq. J3J) (see e.g. ref. |§])- 
Instead of establishing a link between them through a 
suitable choice of the squared matrix element |Mi/| 2 , we 
will try to obtain the formula (TJ3) starting from a suitable 
ansatz which will enable us also to recover quantum statis- 
tics effects (Bose-Einstein or Fermi-Dirac correlations) in 
a natural way. This can be accomplished because a finite 
volume is involved in Eq. (Tj^) unlike in Eq. (|2*|). 

We assume that, as a result of non perturbative QCD- 
driven evolution, the cluster state develops uniform pro- 
jections over the multihadronic Fock space states defined 
by its volume and compatible with its quantum numbers. 
Thus, if |i) is a properly normalized asymptotic state char- 
acterized by the mass, spin and quantum numbers of the 
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cluster and (/| an asymptotic multihadronic final state, 
the rate ry into the final state / is written as: 

^/ = KiW)| 2 (3) 

where W is an effective transition operator proportional 
to the projector over the Hilbert subspace defined by all 
stationary multihadronic states \hy) within the cluster, 
namely: 

W = Y J \hv){h v \f) = Pvn (4) 

hv 

where r) is an operator depending on strong interaction 
symmetry group invariants (Casimir operators) such as 
mass, spin, isospin, charge etc. The state \hy) will be as- 
sumed as a confined stationary free particle state within 
the cluster, with fixed or periodic boundary conditions; 
the inclusion of all resonances as independent states al- 
lows to take into account a part of the interaction between 
strongly stable hadrons and this is the reason of the usual 
expression ideal hadron-resonance gas |15|. 

The operator W is a peculiar one because it is depen- 
dent on the shape and volume of the cluster, which in 
fact pertain to the initial conditions. If cluster's quantum 
numbers coincide with those of the initial colliding system, 
(only one cluster is produced) W should commute with all 
conserved quantities in strong interaction to ensure the 
due selection rules, though this may not be necessary if 
many clusters are produced. 

The commutation requirement is fulfilled for all inter- 
nal symmetries, charge conjugation and for angular mo- 
mentum and parity provided that the cluster has spheri- 
cal shape (see Appendix A). On the other hand, W does 
not commute with energy and momentum as translational 
symmetry is broken by the assumption of a finite vol- 
ume, hence a violation of energy-momentum conservation 
of the order of the inverse of the cluster's linear size is 
implied. However, as it will become clear in the follow- 
ing, momentum-integrated rates in fact get contribution 
only from states fulfilling energy-momentum conservation; 
otherwise stated, finite volume introduces a smearing ef- 
fect on energy and momentum which is washed out af- 
ter kinematical integrations. It should also be pointed out 
that viewing a short-lived object such as a cluster as an 
asymptotic state with definite total energy and momen- 
tum is certainly an approximation and a slight violation 
of energy-momentum conservation is not to be taken as a 
serious awkwardness. Problems may arise only in handling 
single-cluster collision events, where final states must have 
the energy and momentum of the colliding system. 

The Eq. can be written as: 

r s = {f\w\i){i\w^\f) = {f\wp i w^\f) 

= h\ 2 (f\PvP t Pl\f) (5) 

where Pi is the projector over the initial quantum state 
and r\i is such that rj\i) — rii\i). In principle, the pro- 
jection is to be carried out onto a state with definite en- 
ergy, momentum, spin (the Pauli-Lubanski vector), parity, 



C-parity (if the cluster is neutral) and internal charges. 
Hence the most general projector to be considered reads: 

Pi = Pp,j,a, w P x Pi,j 3 Pq (6) 

where P is the four-momentum of the cluster, J the spin, 
A the helicity, tt the parity, x the C-parity, / and I3 the 
isospin and its third component and Q = (Qi, . . . , Qm) 
a set of M abelian (i.e. additive) charges such as baryon 
number, strangeness, electric charge etc. Of course, the 
projection P x makes sense only if J3 = and Q — 0; in 
this case, P x commutes with all other projectors. 

A state with definite four-momentum, spin, helicity 
and parity tranforms according to an irreducible uni- 
tary representation v of the extendend Poincare group 
10(1,3)^, and the projector Pp, can be written by us- 
ing the invariant, suitably normalized, measure fj, as: 

Pp,j,x,« = ±J2 6hxiu f A ^9z) D v \g z )\ U{g z ) (7) 
z =i,n J 

where z is the identity or space inversion g z € 
10(1, 3)^_, D v {g z ) is the matrix of the irreducible repre- 
sentation v the initial state i belongs to, and U(g z ) is the 
unitary representation of g z in the Hilbert space. Simi- 
lar integral expressions can be written for the projectors 
onto internal charges, for the groups SU(2) (isospin) and 
U(l) (for additive charges). Although projection operators 
cannot be rigorously defined for non-compact groups, such 
as Poincare group, we will maintain this naming relaxing 
mathematical rigour. In fact, for non compact-groups, the 
projection operators cannot be properly normalized so as 
to P 2 = P and this is indeed related to the fact that |i) 
has infinite norm. Still, we will not be concerned with such 
drawbacks thereafter, whilst it will be favourable to keep 
the projector formalism. Working in the rest frame of the 
cluster, with P = (M, 0), the matrix element D v ^{g z )\ 
vanishes unless the Lorentz transformations are pure ro- 
tations and this implies the reduction of the integration in 
from IO(l,3) T to the subgroup T(4) <g> SU(2) <g> Z 2 (see 
Appendix B). Altogether, the projector Pp t j t \^ reduces 
to: 

Pp,i,A,, = 7^j4 / d ^ e' P x C/(T(x)) 

x (2 J + 1) J dR D J (R)** U(R) ' + 7r 2 t/(n) (8) 

dR being the invariant SU(2) measure normalized to 1. 
The invariant measure d 4 x of the translation subgroup 
has been normalized with a coefficient l/(27r) 4 in order to 
lead to a Dirac delta, as shown below. This is indeed the 
general expression of the projector defining the proper mi- 
crocanonical ensemble, where all conservation laws related 
to space-time symmetries are fulfilled. 

Hereafter, we will confine ourselves to clusters with 
fixed energy, momentum and abelian charges while conser- 
vation of angular momentum, isospin, parity and C-parity 
will be disregarded. This is expected to be a very good ap- 
proximation in high energy collisions, where many clusters 
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are formed and these latter constraints should not play a 
significant role |15Ul()j . On the other hand, they cannot be 
disregarded in very small hadronizing systems (e.g. pp at 
rest and, in such circumstances, the full projection 
operation in Eq. JSJ should be carried out. As has been 
mentioned in the introduction, the set of states with fixed 
values of energy, momentum and abelian charges is de- 
fined microcanonical ensemble as well and we will stick to 
this convention. 

Dealing with clusters with an unspecified value of an- 
gular momentum isospin and parity means, from a statis- 
tical mechanics point of view, that all possible projections 
over definite values of those quantum numbers occur with 
their statistical weight. In other words, we shall sum over 
alU, A, /, ^3 , 7r, x, which amounts to simply remove the rel- 
evant projection operators in virtue of the completeness 
relations such as, for instance: 



5^(2 J + 1) f dR £> J (R)f U (R) = I 



(9) 



In this case, the projector operator onto the initial state 
reduces to the more familiar form: 



PpP 



1 



(2tt)< 



d xe 



lP-Xg— lPop-X 



^ ) J —TV 



(2 



where — (n, . . . , 7r) and the group generators P op and 
Qo P have been introduced. The appearance of Dirac and 
Kronecker deltas in Eq. (|10|l reflects the abelian nature of 
the leftover space-time translations and U(l) groups. By 
using the latter expression of the projector P^, the Eq. 10} 
and inserting two identity resolutions, Eq. © turns to : 

r f= i^i 2 E E(/im(^i/')(/"i^x^i/) 

h v h' v f'f" 

x(f'\6 4 (P-Po P )8 Q ,Qjf") (11) 

and, taking /', /" states as energy- momentum and charges 
eigenstates: 

r f= E E</IM(MW'l^}(^l/> 

h v h' v /' 

l4/D D - - (12) 



xfi*(P-Pf>)6 QtQfl 



that is: 



r f = h\ 2 E | ^2(f\hv)(hv\f) 5\P - Pf>)S QMf , 

f h v 

(13) 

Multiparticle states in the Fock space are characterized 
by a set of integer occupation numbers for all the species 
and for all the kinematical states. This also applies to 
the general state \hy} as long as it represents, as it has 
been assumed, free hadron and resonance states within 
the cluster, so one can write \hy) = \{Nj}ky) where 



{Nj} = (Ni, . . . , Nk) is a K-uple of integer numbers one 
for each hadron species j and fcy denotes a set of kinemati- 
cal variables, depending on the spacial region with volume 
V, describing the state of the TV = Ni + N 2 + ■ ■ ■ + Nk 
particles. Similarly, we can rewrite the states belonging 
to the complete basis as |/) = \{Nj}k) where now k is 
meant to be a set of proper momenta and polarizations. 
Note that the expression (|13|) allows transitions to states 
|/) with energy-momentum different from P, unless the 
volume is infinitely large. This tells us, as has been men- 
tioned, that the energy-momentum spread is of the order 
of the inverse of the cluster's linear size. 

To further develop Eq. I|13|l we shall assume that: 

({N j }k v \{N j }k)=0 XNj^Ni Vj (14) 

Hence, it is required that states with different particle 
composition, either within the bounded region or in the 
whole space, are orthogonal. Indeed, there are two con- 
traindications to this assumption. The first is of a more 
fundamental character: in relativistic quantum field the- 
ory a condition like (|14|) cannot be exactly true as sta- 
tionary states localized in a finite region are not eigen- 
states of the properly defined particle number operator 
(localization involves the creation of particle-antiparticle 
pairs). However, this effect is relevant if the size of the re- 
gion is lower than the Compton wavelenght of the particle 
1/m, which is at most (for pions) « 1.4 fm, corresponding 
to a volume of « 3 fm 3 ; for all other hadrons, this vol- 
ume is significantly smaller. Henceforth, we will assume 
that volumes to be dealt with are larger (not too much 
though) and will take a non-relativistic quantum mechani- 
cal treatment as a good approximation. The second is con- 
cerned with strongly decaying resonances, which, in prin- 
ciple, should not be orthogonal to the states of their decay 
products; however, we have assumed that resonances are 
to be treated as independent states, so the orthogonality 
relation is correct in the framework of the ideal hadron- 
resonance gas. 

With these two caveats in mind, we proceed to calcu- 
late the total rate of some channel, i.e. a multihadronic 
configuration {Nj}, by summing over the physical observ- 
ables fc, being |/) = \{Nj}k). Applying the sum to the 
right hand side of Eq. I|12|) , taking into account the con- 
dition l|14|) and the completeness of the set \{Nj}k), one 
obtains: 

J2({Nj}k\{Nj}k v ) 

k j k 

x<{Ar i }^|{Ar i }fe) = (h v \ti v )l[5 Nj ^= 5 hvK ]Js Nj ^ 

j j 

(15) 

and, therefore: 

k 

= h\ 2 E E i({^}fcW-}Mi 2 <s 4 OP - Pf) *q,q,,(i6) 

k' kv 
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The above equation (and, maybe more apparently, 
Eq. (|19|l below) shows that only kinematical states ful- 
filling energy-momentum conservation contribute to the 
total rate of a channel even though the transition to final 
states with Pj ^ P is allowed, as discussed. 

We are now in a position to recover an expression 
like 0) mentioned at the beginning of this section. More 
specifically, we can prove that the right hand side in 
Eq. I|16fl is \r]i\ 2 times the usual expression of the proba- 
bility of the multihadronic configuration {Nj} to occur in 
the microcanonical ensemble of an ideal hadron-resonance 
gas with four-momentum P, charges Q and volume V, as 
long as the aforementiond relativistic quantum field effects 
are disregarded. We will show this first in the simple case 
of a channel with all different particles, i.e. Nj < 1 Vj; 
the case of identical particles will be handled in the next 
section. The scalar product in Eq. i|16|) factorizes, so that: 



N 

^ i({^-}fc'i{^-}fc y )i 2 = n e kp»^>i 2 ( i7 ) 

fey i — 1 ki ,Ti 

where i = 1, . . . , N is the single-particle index. The vari- 
able p is a momentum whilst k denotes three variables 
defining the state of the particle within the region with 
volume V (e.g. a plane wave vector for a rectangular box 
or energy and angular momenta for a sphere). The vari- 
ables and r, labels different polarization states of the 
particle and may refer to different projections of the spin 
(or the helicity); we will assume that the transformation 
from t to o is unitary. As long as |k, r) is a complete set 
of one-particle states in the region with volume V, as a 
consequence of the completeness of the states \hy}, it can 
be shown that, in the non-relativistic quantum mechanics 
approximation (see Appendix C): 



EKp*>I 2 = 



V 



k,r 



(2tt) 3 



(18) 



Thus, taking into account that ^2 k , — [Q i=1 Yla I d 3 Pi 
2 and Eq. lfP7|l. Eq. {03) becomes: 



{Nj} 



(27T) 

\r)i\ 2 n {Nj } 



V N r N r i 

' i=i J j t 



(19) 



where charge conservation . Njqj = Q is understood. 
Therefore, the rate P{n 3 } is proportional to the usual ex- 
pression of the phase space volume or density of states per 
four- momentum cell f2{N-} of the multihadronic configu- 
ration {Nj} in the microcanonical ensemble of the ideal 
hadron-resonance gas. It is to be emphasized that this 
formula, which has been used by many authors in the 
framework of the statistical model, is not the most general 

2 The notation [Y[ i J d 3 pi] stands for the integral operator 
J d 3 pi . . . J d 3 pjv and it is understood to act on its right hand 
side argument 



though, as all particles must be different. Therefore, it cor- 
responds to assuming the classical Boltzmann statistics. 
We will see in the next section that, if quantum statistics 
are taken into account, the integral in Eq. (|19|) is indeed 
a single term of an expansion. 

Even though our derivation might look unnecessarily 
elaborated, the TV-body relativistic phase space volume in 
Boltzmann statistics Eq. (|19fl has been recovered start- 
ing from purposely built quantum mechanical transition 
probabilities without invoking any time reversal argument 
or averaging procedures like in previous treatments |9lll(Jj . 
Furthermore, it should be emphasized that this derivation 
is more general than previous treatments because we did 
not consider particle states within the cluster as energy- 
momentum eigenstates. In fact, Eq. Ijl9|l is obtained in a 
traditional approach |1 lj by working out the expression: 



Nj},{Nj} BtB . 



(20) 



with the key assumption that particles within the cluster 
have indeed definite four-momenta and so does the whole 
multiparticle state: 



-Ptatc — ^ 



Pi 



(21) 



and finally replacing the sum over particle states in the 
cluster with an integration: 



E 



V 



(2tt) 3 



d 3 P 



(22) 



Altogether, this approach can be a good approximation 
only for very large volumes, because only for large vol- 
umes can the uncertainty on energy-momentum entailed 
by localization within the cluster be negligible. For smaller 
volumes, the localized multiparticle states, that we have 
denoted as |/iy), are not eigenstates of energy-momentum 
and their spread in energy-momentum cannot be ne- 
glected. It is worth making a rough estimate of how large 
the volume ought to be for the traditional approach to be 
valid. This can done in two ways: requiring that the un- 
certainty in momentum for a single-particle localized state 
is not larger than order of, say, 10% or arguing that the 
approximation (|22fl is indeed a good one provided that 
the number of phase space cells is at least of the order 
of 10-100. Working out (|22|l or using the indeterminacy 
principle, it turns out in both cases that the cluster's lin- 
ear size should be larger than ss (6-10)/p, where p is the 
typical momentum of the particles at hadronization. Since 
this is of the order of some hundreds MeV, the linear size 
must be of the order of, say, 3-10 fm, which is consistently 
larger than the limit set by the aforementioned condition 
on the Compton wavelength of the particles, of the or- 
der of a fraction of fermi. Therefore, the requirement on 
the volume for the validity of the traditional treatment is 
more stringent than that needed for the present one. 

Although in the case of Boltzmann statistics, the tra- 
ditional and the present approach lead to the same ex- 
pression Eq. I|19[) for the A^-body relativistic phase space 
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volume, different expressions are found in the case of quan- 
tum statistics, that is with identical particles, as it will be 
shown in the next section. 



3 Identical particles and cluster 
decomposition 

If there are identical particles in the channel, the equa- 
tion l(X6"|l holds but Eq. (fPTjl does not and has to be mod- 
ified. For sake of simplicity, we will start with the case of 
only one kind of particle in the channel and assume that 
charge conservation is fulfilled. As has already mentioned, 
relativistic quantum field effects will be disregarded, nam- 
ley cluster's size is assumed to be significantly larger than 
the Compton wavelenght of the particle. The correspon- 
dance between Fock space and multiparticle tensor space 
requires the identification: 



where, in the last equality, we have redefined the dummy 
permutation indices p, q as their inverse and multiplied 
each term by a factor n\\ . . . Um-/N\ in order to avoid 
multiple counting of the symmetric (antisymmetric) basis 
tensors \hy) when the sum over all possible vectors k and 
polarizations r is carried out. Finally, taking into account 
that also and Ti are dummy indices, one sum over per- 
mutations can be trivially performed and we are left with 
the transformation: 

E E i mmN^kv) i 2 -> i e 

k k v r 
N . 

X [ll / ^P 1 ] El (Pi Iki^i) < k i T i I Pr(2) Cr(i) ) (27) 



\{Nj}k v ) 



being r = p x q and x(0 = x(P = xip l )x{o) = 
x(p)x{q)- The inner sums in the above equality yield (see 
Appendix C): 



p v 



1 



|l%)(l)Tp(l), . . . , kp(N)T p (N)) (23) 



k;T, 



where p is a permutation of the integers 1, . . . ,N and x(p) 
its parity; the rij's are the number of times a given vec- 
tor ki recurs in the state with X)f=i ni = N; b — for 
bosons and b = 1 for fermions. As there is only one parti- 



(28) 

so the following expression of the phase space volume Qn 
for N identical particles is obtained: 



cle species, the phase space volume ^j^} can be denoted Q N — 'xfj) /" d 3 pi . . . d 3 pAr £ 4 (P — P/) 

with Hat and can be calculated bv using En. fTTHl. Renlac- > iV! „ „.J 



with Qm and can be calculated by using Eq. Ijl6|l . Replac- 
ing \{Nj}k') with: 

\{Nj}k) -> E ^tfT IPpW^pU)' ■ ■ ■ :Pp(JV) cr p(JV)> ( 24 ) 



p 



TV! 



i r 

(2 7 r)3i y dXC 



(29) 



similarly to Eq. and dividing by 1/JV! in order to 

avoid multiple counting of (anti-) symmetric basis tensors 
when integrating over all possible momenta, we find: 



r, 



JV 



N 

[nE 



6 4 (P-P f 



El 



x(p) 1 



1 



VJY!ni!...nji/! 



(Pi^i, • 



tp(l)T p (l), 



(25) 

In Eq. 1251) and hereafter Pf must be understood as the 
sum of the four-momenta of all particles in the channel. 
The last factor in the above equation can be worked out 
as follows: 

E I E ^ ^J,,,^, ^. ■ ■ • iWpm. ■ 

fey p 

x 1 



^ N\n x \ 

fey 
X(pi0"i, 



..n M 

k p(i) r p(i) 

N 



iPi^i, 



)( k 9(l) T 9(l) 



p.q 



i—1 kiTi 



Hence, the phase space volume of N identical particles 
is given by the sum of N\ terms and it is thus enhanced 
or suppressed with respect to the case of distinguishable 
particles. As it will be proved in the following, this ef- 
fect is owing to the finite volume and, thereby, this model 
naturally accounts for Bose-Einstein and Fermi-Dirac cor- 
relations. 

To develop Eq. (|29[) . it is useful to recall that any per- 
mutation r of N integers can be uniquely decomposed into 
the product of cyclic permutations, that is r = c± ... ch ■ 
Let n be the number of integers in each cyclic permutation 
and let h n be the number of cyclic permutations with n 
elements in r so that JDnLi = N. The set of integers 
h%, . . . , hjy = {h n }, with J2n°=i = H, is usually defined 
as a partition and different permutations having the same 
structure of cyclic decomposition, that is the same num- 
ber of integers for each q (i.e. the same partition), belong 
to the same conjugacy class of the permutation group Sat. 
The crucial observation is that each term in Eq. (|29[1 is in- 
variant over a conjugacy class, or, in other words, depends 
only on the partition {h n }; this happens because different 
permutations in the same conjugacy class differ only by 
a redefinition of the integers 1, . . . , N and this is just a 
change of the name of the dummy integration variables 
and sum indices in Eq. I|29|) . The number of permutations 
of Sjv belonging to a given conjugacy class is a well known 



(26) number [TZ], namely Nl/Jj 



N 



l h n \. Furthermore, for 
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the representative cyclic permutation c with n elements 
(l,...,n): 



e n^-^+D 



(30) 



so that: 



N 



H 



e n^,=n e n<^>=(^+i) H 



r(0 ±± ±± 

fTi,...,crjv i=l l^l(Ji,...,(T n — 1 



where r = c\ . . . Ch and ni is the number of integers in the 
cyclic permutation q. By defining: 



F„ 



n 



i 



(27T) 



d 3 xe ix-(p Cl(ll) ~ Pl 



(32) 



and taking into account that x( c = ( — 1)"' +1 , we can 
finally rewrite Eq. I|29ll as: 

H N = ^ ^N({h n }) = 

{hn} 

( T 1) N + H (2J+1) H 



E 

{h n } 



N 

[n 



H 



d 3 Pj 



S\P-P f )]jF ni 



(33) 



where the upper sign applies to fermions, the lower to 
bosons. Therefore, the phase space volume for a channel 
with N identical particles consists of a large number of 
terms /2jv({ft. n }), each corresponding to a partition {h n }, 
what is usually called in statistical mechanics cluster de- 
composition 3 . 

In the large volume limit, the dominant term is the 
one with the highest power of V and this corresponds 
to the partition (h\, h 2 , ■ ■ ■ , hx) — (N, 0, . . . , 0), i.e. the 
identical permutation. In this case there are N factors 
F\ = V/(2tt) z and the whole term reads: 



n N (N,o,...,o) 



iV(2J + 1)- 


n i r 


i (2tt) 3 . 


mj 



where Pf = 2p2 + P3 
gration variables p' = 



- . . .pN- Introducing the new inte- 
2p 2 the energy term 2e 2 becomes 



\/p' 2 + (2m) 2 and Eq. I|35l) can be rewritten as 
/2jv(iV-2,l,0,...,0) 



(Tl) r V(2J+l) i n-i 
2 4 (iV - 2)! L (2tt) 3 . 

r - 

x I d 3 V 'd 3 V3 ...d 3 VN 6 4 (P-p'-J2Pi) (36) 

i=3 

Aside from the sign and an overall normalization factor 
1/16, this term corresponds to the Boltzmann limit 1)35(1 of 
the phase space volume for a set of N— 2 identical particles 
plus a new particle (labelled with a prime) obtained by 
clumping particles 1 and 2 into a lump with a mass twice 
the mass of 1 and 2 and the same spin. 

Actually, this kind of interpretation holds for all of 
the terms in Eq. (|33|l . In fact, in the large volume limit, 
each F n implies the elimination of n — 1 integration vari- 
ables through the appearance of Dirac deltas, while a sin- 
gle V/(2tt) 3 factor is left because of the cyclic structure of 
the permutation, namely: 



p 

± 71 1 



V 



"1-1 



(2tt) 3 J-l 

V ' ii=l 



n ^(p^ -p*i+i) 



(37) 



Then, after trivial integrations are carried out in Eq. I|33|) . 
the Dirac delta forcing conservation of four-momentum 
turns into 5 A (Pi - n x p x - n 2 p ni +i ■ ■ ■ - n H p nH -i+i) an d 
new integrations variables can be introduced: 



pi = raipi p 2 



«2Pm+l 
as well as new energies: 



p'h = n H Pn H -i+i (38) 



nie = ni\/p 2 + m 2 = \J p' 2 + (nimy 



(39) 



Therefore, the term corresponding to the partition {h n } 
can be written as: 



d 3 pi 



.dVsV-Pf) (34) Q N {{h nj }) 



(Tl) 



N+H 



V(2J+ 1) 



(2tt) 3 



H 



which is the phase space volume for a set of N iden- 
tical particles in the classical Boltzmann statistics; we 
have indeed recovered the phase space volume quoted in 
Eq. Q19f) . All other terms of the expansion in Eq. (|33|l have 
a lower power of V. The next-to-leading term correspond 
to the conjugacy class of permutations with one exchange 
and N — 2 unchanged integers, i.e. (hi, h 2 , /13, . . . , ft. at) = 
(N — 2, 1, 0, . . . , 0) In the large volume limit it is easily 
seen, looking at Eq. $2$, that F 2 -> <5 3 (p 2 - p 1 )V/(2w) 3 
and the whole term thus reads: 



/d3pi 



a 6 p' H s*(p 



X>0 

1=1 



(40) 



where particles are now clumped into H lumps with mass 



equal to nim and spin J. Since 0^=i 
above equation can be written also as: 

(=fl) N + H rV(2J- 



,3 _ 



n 



N 

71=1 1 



,3h„ 



the 



n N ({h nj }) 



(2tt) 3 



H 



(Tl) \V(2J+l) 



x 



2(JV-2)! 
Jd 3 p 2 ...d 3 p N 5 \P - P f ) 



(2tt) 3 



N-l 



x |dy 1 ...dy ff <s 4 (p-X>0 (4i) 



/=i 



(35) 



3 The term cluster in this context has nothing to do with our 
previous definition of an individual hadronizing source 



We can thus conclude that the general term relevant to 
the cluster decomposition of the phase space volume of a 
set of N identical particles can be obtained by calculating 
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the phase space volume, in the Boltzmann statistics, of 
a suitable set of lumps having as mass multiple integer 
values of m and spin J, weighted by an overall coefficient 
of ( T 1) N ~ H / H n n 4h " . Note that the factors l/h n \ already 
take into account the identity of the lumps. 

After having inferred the expressions of the phase 
space volume of a channel with TV identical particles, the 
generalization to a channel {Nj} (see Sect. 2) with an ar- 
bitrary number of groups of identical particles for each 
species j is rather straightforward and can be achieved 
by going along the previous arguments. Thereby, the fol- 
lowing equations are obtained which are extensions of 
Eqs. 123), 1231 respectively: 



k ky 



Ni. 

3 J rjES Nj 



X 



n {Nj} = E / d3 Pi ■ ■ ■ d3 P^ s4 ( p - p /) 



n e 



NA 



b N i s 3 r 

TT -jttt^ / d 3 xe IX - (p ^^'^ ) (43) 



and 



Q 



d 3 Pl ---d 3 pjv 5\P-P S ) 



n e ™!TZ:r n <«> 



' ' TT J, i 

3 {h nj } 11„ 3 =1 



with Hj = J2n-=i h n:i an d Nj = Yln-=i n jh nj ■ The above 
expression is the most general for the microcanonical 
phase space volume of the multihadronic channel {Nj} in 
the ideal hadron-resonance gas framework with full quan- 
tum statistics and generalizes the expression obtained in 
ref. [HI: 



{Nj} 



TT V (Tl) Ni+H > - 

lln 3 =l 11 j """■j ■ 



3 {h nj } 
H 



[n 

U=i 



V{2Jj + 1) 



(2tt) s 



d 3 P ; 



Hi 



<5 4 (P i -E^)(45) 



where, for a set of partitions . . . , {h nK }, the four- 

momenta are those of lumps of particles of the same 
species j (Hj in number) with mass rijmj and spin Jj. 
While Eq. H44fl is indeed the correct expression of the mi- 
crocanonical relativistic phase space volume of the chan- 
nel {Nj} in the statistical hadronization model, Eq. 
turns out to be a special case of Eq. (|44l) and in fact can 
be derived from it by replacing F ni with their limiting 



expressions Ij37(l . Thus, the expression Q45J1 is a good ap- 
proximation of 14411 only in the limit of large volumes. In 
fact, the derivation of the relativistic phase space volume 
in Eq. (|45|) in ref. mi was based on the assumption that 
particle states within the cluster are energy-momentum 
eigenstates, which is a good approximation only for large 
volumes, as extensively discussed at the end of the previ- 
ous section. For the expression Eq. (|44|l to be valid, we now 
just need the linear size of the cluster must be sufficiently 
larger than the Compton wavelenght of the involved par- 
ticles in order to neglect relativistic quantum field effects. 

The fact that the leading Boltzmann terms in the gen- 
eral expression (|44Jl and the approximate one Q45J1 are the 
same, as already pointed out in Sect. 2, reduces the ac- 
tual numerical impact of this generalization on many, yet 
not all, observables. For instance, at the actual tempera- 
ture values of about 160 MeV found in previous analyses 
of many high energy collisions in the canonical ensemble 
quantum statistics corrections on average particle 
multiplicities turned out to be significant for pions only 
(about 10%), whilst they can be neglected for all other 
hadrons. Therefore, as long as average multiplicities are 
concerned, the calculation can be done within Boltzmann 
statistics and the difference between the correct formula 
and the approximate one is almost irrelevant. Thus, the 
only effective requirement on cluster size for the validity 
of all performed analyses in high energy collisions P]E] is 
that it must be larger than Compton wavelenght of parti- 
cles (at most 1.4 fm) and this is always met. 

Even though average multiplicities are essentially un- 
affected in most practical cases, there are other observ- 
ables which are sensitive to quantum statistics effects and 
for which the fully correct calculation of phase space vol- 
ume^] is compelling, e.g. Bose-Einstein correlation spec- 
tra and multi-pion exclusive channel rates. 



4 Microcanonical partition function 

The overall phase space volume of the ideal hadron- 
resonance gas is obtained by summing &{n } over all al- 
lowed channels: 



n = E Q {Nj} S Q,Q {Nj} 

{N 3 } 



(46) 



As f?{Nj} — r{Nj}/\r]i\ 2 , fi can be expressed on the basis 
of Eq. Hll|) after having removed the two identity resolu- 
tions in /' and /": 



= E E (/IM(W|<5 4 (P-Pop)^Q,Q OP l^)(^l/> 
/ hv.h' 



E(VI<5 4 (P-Pop)^Q, Qo pl^ 



(47) 



The last expression makes it apparent that the definition 
of fl as the microcanonical partition junction is an appro- 
priate one. If the sums in Eqs. (|46|l and (|47[) are performed 
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over all channels regardless of their charge, the obtained 
quantity is defined as grand-microcanonical partition func- 
tion: 



(48) 



Throughout this section we will confine ourselves to the 
latter, rather than the to properly defined microcanonical 
partition function, in order not to bring along a cumber- 
some formalism. This limitation shall not affect the gen- 
erality of the expounded arguments and the extension to 
the case of constrained charges is indeed straightforward. 

We have seen in the previous section that Eq. 14411 
is a correct generalization of Eq. (|45|l for finite volumes. 
Likewise, Eq. (|48|l is a generalization of the expression 
quoted in previous literature [lllTfj : 



states 



(49) 



In fact, Eq. I|49|) is a straightforward consequence of 
Eq. (|48|l if \hv) is an eigenstate of energy-momentum. 
However, we have already emphasized that \hy) is a local- 
ized state and its four-momentum is a well defined quan- 
tity only in the large volume limit, as has been discussed at 
the end of Sect. 2. Thus, Eq. l|4*9"|l is consistent only if the 
cluster is sufficiently large, whereas the Eq. I|48[l is always 
a well defined one. On the other hand, a closed analytical 
integral expression for the correct (grand-)microcanonical 
partition function cannot be written. The best one can 
do is to decompose it as a sum over channels, as in 
Eq. H4fi|) . and calculate numerically the Q{n -}' s according 
to Eq. I|44|) . which is a formidable task indeed. Conversely, 
Eq. (|4*9l does lead to a closed integral expression, which 
can be obtained by firstly expanding the Dirac delta in 
Eq. (|49J) as a Fourier integral: 



6 4 (P-Pst,te) 



1 



{2-kY 

and reexpressing Eq. ll4^|) as 



d 4 y e 1 



(-P-Pstatc)-y 



(2*) 



-inj h pj h -y 



(50) 



(51) 



{n jh } j,h 



The sum over states is in fact a sum over all possible 
occupation numbers of each phase space cell. The calcula- 
tion now proceeds by taking advantage of the commutabil- 
ity between sum and product in l|51|l . However, unlike for 
fcrmions for which rijh = 0, 1 only, the sum over occupa- 
tion numbers does not converge to a finite value for bosons 
as Ujh runs from to oo. The convergence is recovered if 
the time component of y is provided with a small nega- 
tive imaginary part — ie. If we introduce such a term in 
Eq. (|51|l the sums can be performed and the result is: 



Q = lim \ dy a / d 3 y e lP y 

e^O (27T)4 J^J J 

xexp [^log(l±e~ ip ^" y ) ±1 

j,h 



(52) 



where the upper sign applies to fermions, the lower to 
bosons. The integrand function is in fact singular for y = 
and the shift of the integration contour in the complex 
plane provides a regularization prescription. The sum over 
phase space cells can be replaced, in the large volume 
limit, by an integration according to Eq. I|22|) . so that i? 
reads: 

-i r+oc— ie p 

fl = lim — — I dy° / d 3 y e lP y 



+o (2tt) 4 



— oo — IS 



x exp 



3 



(2tt)3 



d 3 p log(l ± e^-yf 1 } (53) 



We will prove in the remainder of this section that the 
closed expression for the grand-microcanonical partition 
function, Eq. i|53|) , can be recovered without invoking (|49H 
and H22f) . but starting from the general expression l|44|) in 
at least two cases: 

1. for Boltzmann statistics; 

2. in a full quantum statistics treatment, by enforcing the 
approximation Eq. I|37|) . namely: 



(2 



^ / dWp-p'^-^p-p') 

V Jv 



(54) 



Henceforth, we will adopt the following shorthand: 

r +oo — is 

dy° I d 3 y (55) 



dV 



1. Let us start by showing that for Boltzmann statistics. 
We have seen in the previous section that confining to 
classical statistics amounts to retain only the first term 
{h nj } — {Nj, 0, . . . , 0) in the general cluster decomposi- 
tion Eq. igD, hence h x = H 3 = Nj, Fx = V/{2n) 3 and 



/~}Boltz 



ff= |d 3 Pl ...dW 4 (p 

1 \V{2J+l 



N 

E 



Pi) 



n 



Njll (2nf 



A', 



(56) 



The Dirac delta in the above equation can be Fourier ex- 
panded, thus, after regularization: 

1 /-+00 — IE p 

*™'> = (w/-„J> e '''/^'-' dV 



JY 



1 \V{2J+l)iN 

3 



«p[-i£»-dll]^ 

1=1 j J y 

-V(2J+l) 



(2nY 



+oo — ie -i 
-oo-ie ,• JV J- 



(2tt) 3 



d 3 p e" 



(57) 

Summing over all channels yields the grand- 
microcanonical partition function: 



fl Boltz = 



(2tt)4 
x exp |^ 



-oo — ie 



d 4 y e 1 



OO — 16 



j 



F(2J ^ + 1) / d 3 p o 



(58) 
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which can be obtained indeed from Eq. I|53|l with the 
boltzmannian approximation: 



log(l±e- ip - 2 ') ±1 ~ e - ipy 



(59) 



This proves the first part of our argument. 

2. If quantum statistics is included, we make the sup- 
plementary assumption, as has been mentioned, that ap- 
proximations H37|) a Pply and, thus, Eq. 144|) turns to 
Eq. |@SJ. Let us first restore pi j — pj./n^ (see Eq. (J3SJ) 
as integration variables and rewrite Eq. (|45|l by plugging 
in the Fourier expansion of the Dirac delta: 



n 



{N 3 -} 



n e c*i) 



3 {h nj } 



H 



n 

i J= i 
i 

Hi 



(2tt) e 



n ni h ' 



+ 00 — IE 

d 4 v e lP -" 



(27T) 
^+oc— is 



3 '{"} II,) : &n ! 

3 lj=l 



— — j / d 4 ye^TT V (Tl) 

" ; 7 x j :/•• :■ II, i ".. I'-.- 



(27T) 



x j| [ g&± ^ / <l ! i»/ ^ 

(,'=1 



(2tt) e 



(60) 



In order to simplify the notation, we introduce the quan- 
tities: 



^ E ^ (y)= Tr / d3pe ' 1,,M 

so that Eq. (|6U[1 can be further written as: 
1 



(61) 



(27T) 



-oo — ie 

d 4 y e lP y 



— oo — is 



i 



(2tt) 4 



+oo — le 



dv^nE n 

j {/i n ,}™j =i 



(T1) K+i)h„ jz 



J 



(62) 



where, in the last passage, we have taken advantage of the 
fact that Zj(n) is constant over a conjugacy class. Also 
note that we have released the upper limit in the sum be- 
cause of the constraint ^ n rij h n . = Nj which effectively 
sets it to Nj. At this stage, the key observation is that we 
can implement this constraint through an integration in 



the complex plane for each species j and then perform an 
unconstrained sum over all h r , . : 



E - E s Sn _ njhn . iNj 

{h„.} fii,...,fioo 



1 

2tti 



dw 



oo oo 



ill 



En,- h n 
w 



n,-=l hn-=0 



(63) 

so that the part of integrand in Eq. 162(1 following the 
j-product sign can be written as: 



1 

2th 



dwj 



OO OO 



3rn E 



t-^-i \(n4-\-l)h n ■ n j n-ihn- 



l 

2th 



y j i»j=i/i 
dw 



■it'- 



ll exp 



"J n j= l 
1 d N > 

F7 ex P 



3m ) 

rij ni h nj \ 

(zpl)^+l 



n-j —1 



rij + 1 



(64) 



By using the explicit expression of Z(„) in Eq. 1|) . the 
series in the exponential of Eq. I|64|l can be summed up 
and this yields, for the phase space volume fi^ N .y. 



Q 



d N > 



(27T) 



-1 /-l-OO— IE 1 



dw 



N, 



exp 



(2J,- + l)V 



(2nf 



d 3 p log(l ± Wje- ip - y ) ±:L (65) 



U! 7 — 



We are now in a position to calculate fl by summing over 
all Nj Vj, according to Eq. I)48|l. The sum over each iVj = 
0, . . . , oo can be performed independently and, noticing 
that each term is the N^ h one of the Taylor expansion of 
the exponential function evaluated at Wj = 1, one obtains: 

1 



a = 



+oo— ie 

4„, n \P-y 



(2tt)< 
x ]^[exp 



d 4 y e 1 

—W 

{2Jj + 1)V 



(27T) 



d 3 p log(l±e- ip2 ') ±1 (66) 



which coincides with Eq. (|53() : this proves our second 
statement. 

The recovery of the known expression of the micro- 
canonical partition function in the two considered cases 
is not surprising as the same holds for the single channel 
phase space volume }. We have seen this in Sect. 2 
where it has been emphasized that I?{j\r-} m Boltzmann 
statistics does not differ from its approximation in the 
large volume limit; and in Sect. 3, where we have seen that 
the £2{Nj} m full quantum statistics l|45|l deduced from 
Eq. (|44|) by enforcing the approximation J37J), was ob- 
tained in the traditional approach JTTJ using Eqs. H2U[) . <|21ll 
and iJUJ. 



5 From microcanonical to canonical ensemble 

What has been done for the grand-microcanonical en- 
semble can be straightforwardly extended to the properly 



F. Becattini, L. Ferroni: Statistical hadronization and hadronic microcanonical ensemble I 



11 



called microcanonical ensemble by adding the further con- 
straint of M abelian charges conservation, like in Eq. ((46(1 
The Kronecker delta can be Fourier expanded: 



M 



XQr, 



(2tt 



where the vector notation <fi = (<^x, • • • ,4>m) has been 
introduced. The reasoning in the previous section, from 
Eq. 1(56(1 onwards, can be easily repeated with the addi- 
tional charge constraint 1(67(1 . under the same conditions 
for the validity of the needed approximation 1(54(1 . One can 
thus arrive at the following expression of the microcanon- 
ical partition function: 



We assume that the above equation has one real solution 
(note that Z(z) is real for real argument, see Eq. J7TJ). 
This must be a timelike four-vector for the momentum in- 
tegration in Eq. 1)70(1 to converge. Therefore, we can set 
(3 — (1/T)u where a is a unit timelike vector and T > 
is defined as temperature, while (3 is usually called tem- 
perature four- vector. It is not difficult to verify that if the 
<jM ^ e i(Q-Q{w 3 })-0 (gy-j cluster's rest frame is chosen, where P = (M, 0), (3 has 

vanishing spacial components and the usual expression of 
the canonical partition function is recovered: 



1 



x exp 



(2tt) 4 + m 
'^(2J j + l)V 



d 4 y e 1 



(2tt) 3 



d 3 p log(l±e~ ip - y - iq -'^) ±1 (68) 



where = (qji, ■ ■ ■ , Oju) are the abelian charges of the 
j th hadron species. Let us perform a rotation in the four- 
dimcnsional complex hyperplane by setting z ~ ly and 
rewrite Eq. (|68|l as: 



f2 = lim 
where 



(2tti) 4 



-\-oo-\e 

d 4 z exp[P • z + log Z(z, Q)] (69) 

oo — ie 



Z(z,Q) = 

[? i2;r; :; 



x exp 



1 

'2n) M 
(2Jj +1)V 



d 3 P log(l ± e -*-p-»h-<t>)±i] (70) 



In Eq. (|70|l it is recognizable the expression of the canon- 
ical partition function 0(7] calculated for a complex four- 
temperature z. The same expression can be obtained start- 
ing from the definition: 



Z(z,Q) = ]T(Me- z - p -<5 Q , Qop |/H 



(71) 



and proceeding in the very same way as for the micro- 
canonical partition function. Particularly, the approxima- 
tions l|54"|l are needed to get to Eq. lj7U|l . 

If the volume and the mass of the cluster are large, 
one can make an approximate calculation of the integral 
in Eq. I|69|) through the saddle-point expansion. The large- 
valued parameter can be either volume or mass provided 
that density M/V is a finite value, which is indeed the case 
of interest in the framework of the statistical hadroniza- 
tion model. The saddle-point four-vector (3 is determined 
by enforcing the vanishing of integrand logarithmic deriva- 
tive for each component /x: 



d 
dzf+ 



[P-z + logZ(z,Q)] 



P, 



z={3 



d 



logZ(/3,Q) =0 
(72) 



Z(Q) = X>v|e 



-H, 



Q,Qo P 



(73) 



where H op is the hamiltonian. Retaining only the lead- 
ing term of the asymptotic expansion, the microcanonical 
partition function can be approximated as: 



/2~exp[P-/3 + logZ(/3,Q)] 1 



1 



(27r) 4 detH(/3,Q) 



(74) 



where H is the Hessian matrix d 2 log Z / dz tl dz v . In the 
cluster's rest frame (3 — ((3, 0) as already pointed out, 
thus, according to Eq. (7UJ|, the derivative 9 log Z/dz l 
with respect to the spacial components of z vanish be- 
cause of odd-symmetric momentum integrands and, con- 
sequently, the Hessian determinant in Eq. I|74|) simply be- 
comes d 2 log Z/d(3 2 = C V T 2 . Altogether, if V is large, the 
microcanonical partition function fl is proportional to the 
canonical partition function Z and we can write: 



n(p, Q) oc 



^• p Z(/3,Q) 



with (3 given by Eq. JEJ, and: 



Z(fl,Q) = J d 4 P9(P°)e- f3P n(P,Q) 



(75) 



(76) 



This equation is indeed an exact one, as can be realized 
from Eq. 1(68(1: the canonical partition function is in fact 
the Laplace transform of the microcanonical one. 

The question arises whether and in which range of val- 
ues of cluster's volume and mass the approximation J73J), 
i.e. the use of the canonical ensemble, employed in sev- 
eral analyses of multiplicities in elementary collisions, is 
a good one for the calculation of relevant physical quan- 
tities. This issue can be tackled only numerically for the 
particular system of the ideal hadron-resonance gas, com- 
paring the exact with the approximate calculation; as has 
been mentioned in the Introduction, this will be the main 
subject of the second paper [T^). 

The way temperature has been introduced starting 
from the microcanonical ensemble in Eq. ((72(1 is rather un- 
usual and deserves some discussion. Through the saddle- 
point relation l(72() , we have defined a temperature by en- 
forcing the known values of energy and momentum of the 
cluster to be what it can be easily recognized as the av- 
erage energy and momentum in the canonical ensemble, 
that is, in the cluster's rest frame where P — (M, 0) and 
0=G9,O): 

M=— logZ09,Q) (77) 
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with T = 1/(3. On the other hand, it is also possible 18 
to extend the relation: 



i _ as 

T~ dM 



— = — — with S = log fl 



(78) 



to the microcanonical regime. This definition gives rise to 
the following equation, by using Eq. (55) l: 



- = -\ 1 

T Q (2tti) 4 



+oo— ie 

d 4 

-oo— ie 



:: z°e Mz °Z(z,Q) 



(79) 



which implies a different definition of temperature with 
respect to Eq. f77)l. At the leading order of the asymp- 
totic expansion of the above integral and the previous 
equation reads: 



I /?'exp[M+^logZ(/3',Q)](Gy(/3)/? 2 ) 
T 



2\-l/2 



exp[M + ^ \ogZ(P, Q)](Cv(^)/3' 2 )- 1/2 



where 8 is the solution of Eq. (|77|) and j3' that of 



(80) 



(81) 



If the system is very large, i.e. in the thermodynamical 
limit, the temperature 1/(3' is much less than M so, ac- 
cording to Eq. IP?!)! 8' ~ 8. 1/T ~ 8' ~ 8 and the two 
definitions coincide, as expected. 

It is worth pointing out that, even in the canonical en- 
semble, for finite volumes, \ogZ is not a linear function 
of V (see Eq. JZOJO and this has the remarkable conse- 
quence that T, in both definitions, is not a function of 
M/V but of M and V separately. Otherwise stated, if T 
and V are used as independent thermodynamical param- 
eters, the mean energy is not an extensive variable as it 
does not scale linearly with V. Of course, this mostly un- 
familiar feature disappears in the thermodynamic limit. 



6 Physical observables 

The comparison of the model predictions with experi- 
mental measurements involves the calculation of quanti- 
ties which can be always written as averages or expec- 
tation values of some operator. For instance, the aver- 
age multiplicity of the j th hadron species in the grand- 
microcanonical ensemble can be written as: 



(N 3 ) 



n 



(82) 



the correlations between j th and fc th hadron species as the 
expectation value of (Nj — (Nj))(Nk — (Nj)) and the prob- 
ability of a single configuration {Nj} as the expectation 
value of S ] y i Ni , . . . , 5ft K N . The analytical expressions of 
sums like that in Eq. (|82[1 can be obtained by multiplying 
^{Nj} by a factor (a fictitious fugacity) Xj powered to iVj 
and taking the derivative with respect to Xj for Xj = 1. 



Therefore, for the average multiplicity of the j th hadron 
species: 



-fi ox- log 2- A ^ Q ^ 



{N 3 } 



The sum on the right hand side can be generalized to all 
species: 



G(A 1) ...,A K )= }Z [] {n ] }II X 

{N 3 } 3 



(84) 



and G can be properly defined as the generating function 
of the multiparticle multiplicity distribution. Note that 
G(l) = n. 

The main advantage of this method of expressing ex- 
pectation values is that the generating function can be 
calculated analytically. By using the expression of ^{n,} 
in Eq. i|tj5fl , the right hand side of Eq. (|84|l can be turned 
into: 



G(Ai, . . . , Xk) 
x ex P 



c + OO — 1£ 



(27T) 

(2tt) 3 



1") J-oo-xe 

d 3 p log(l ± Xje^P-y f 1 } (85) 



and similarly in the canonical case. Now the expectation 
value of any operator can be calculated from the generat- 
ing function by applying many times the differential op- 
erators Dj = Xjd/dXj. In fact, according to Eq. JHIJ, for 
the M th power of Nj 



3 ' 

M 



(N 



sin*. 



i=l 



d 

dX~ 



G(Ai, . . . , Xk) 



x=i 



(86) 



Then, since the operators Dj and commute, we can 
write formally, for any function of N\ , . . . , Nk ■ 



1 



(F({Nj})) = -F(D 



,D K )G{X, 



A=l 



(87) 



Thereby, analytical expressions of various observables can 
be inferred. For instance, by using Eq. (|86|) with G given 
by H85J) . the average multiplicity of the j th hadron, in the 
Boltzmann statistics limit, turns out to be: 



1 (2Jj + 1)V 
Q{P) (2tt) 3 



d 3 p n(P-pj) 



where fi(P) is given by Eq. (55} • It is worth remarking 
that, since Q(P — pj) vanishes when (P — pj) 2 < 0, the 
integration in momentum is cut off when, in the cluster's 
rest frame, the energy of the particle exceeds the clus- 
ter's mass, as it should naturally occur in a microcanonical 
framework. 

Despite their simple appearance, expressions like (|88J) 
are extremely hard to calculate analytically. In fact, the 
whole issue of providing closed formulae of multiplicities, 
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correlations etc. reduces to the calculation of the gener- 
ating function in Eq. 1)85 jl. However, an explicit solution 
of that four-dimensional integral is known only in the two 
limiting cases of ultrarelativistic (vanishing masses) and 
non-relativistic gas JU| • F° r the relativistic gas with mas- 
sive particles, which pertains to the hadronic system, no 
closed formula useful for numerical evaluation has ever 
been obtained, not even as a series. Therefore, the only 
practicable way of calculating averages within the micro- 
canonical ensemble is to evaluate ^{Nj} integral expres- 
sions like (|45J) (which is in turn made up of integral terms 
like (jnnjl) an d sum over all possible channels. However, 
also those integrals have been solved analytically only in 
the aforementioned two limiting cases because the func- 
tions to be dealt with are essentially the same. Several 
authors have tried approximations |2()j but in most cases 
it is difficult to keep the error under control so that, at 
some fixed order truncation of the expansions, the rel- 
ative accuracy may vary from some percent to a factor 
of 10 Thus, the problem of exploring hadronic mi- 
crocanonical ensemble can be attacked only numerically 
through Monte-Carlo integration. This has been done by 
Werner and Aichelin in a quite recent paper with a method 
based on the Metropolis algorithm 0. In the next paper, 
we will present a full numerical calculation for the ideal 
hadron-resonance gas which exploits a modification of that 
method, very effective for large clusters, taking advantage 
of the grand-canonical limit of the multiplicity distribu- 
tions as proposal matrix in the Metropolis algorithm. 

7 Summary and outlook 

This paper is the first of a series of two devoted to the 
study of the microcanonical ensemble of the hadron gas, 
which is the most fundamental framework for the statis- 
tical hadronization model. In this work we have mainly 
developed the analytical formalism, while numerical cal- 
culations will be the main subject of the second paper. 
The main achievements can be summarized as follows: 

1. We have provided a consistent formulation of the sta- 
tistical hadronization model starting from purposely 
defined quantum transition probabilities. This formu- 
lation is much easier to handle than previous ones 
based on time-reversal arguments and ^-matrix ele- 
ments averaging and allows to calculate any final-state 
observables more straightforwardly. Furthermore, it is 
easier to extend it to the case of angular momentum 
and parity conservation, whenever needed. We think 
that this formulation clarifies once more that it is pos- 
sible to account for the observed statistical equilibrium 
of the final state hadronic multiplicities as a result of 
prehadronic cluster decays, without invoking a ther- 
malization process driven by collisions between formed 
hadrons. 

2. We have worked out the rates of exclusive channels 
neglecting angular momentum, parity, isospin and C- 
parity conservation (which are important only for very 
small hadronizing systems) and recovered known ex- 
pressions in the statistical model. We have obtained 



an expression of the phase space volume in full quan- 
tum statistics as a cluster decomposition, Eq. 1)44(1 . 
generalizing previous ones |llj which are valid only 
asymptotically, i.e. in the limit of a very large cluster 
(in practice with a linear size roughly larger than 3-10 
fm). This expression is valid provided that relativistic 
quantum field effects are neglected, i.e. the hadroniz- 
ing cluster should be sufficiently larger than Compton 
wavelenghts of the hadrons. 
3. We have shown analytically how the canonical ensem- 
ble can be obtained as an approximation of the micro- 
canonical ensemble for large volumes and mass of the 
cluster. 

In the second forthcoming paper |12| . the numerical in- 
tegration of the microcanonical expressions obtained in 
this paper will be carried out by means of a Monte-Carlo 
method. This will enable a detailed comparison with the 
canonical ensemble and to establish the range of validity of 
the latter, which has been used in the actual comparisons 
with measured hadronic multiplicities PEj. Besides, the 
implementation of a reasonably fast and reliable Monte- 
Carlo algorithm for the microcanonical hadronization of 
single clusters in high energy collisions is a decisive step 
for further tests of the statistical hadronization model. 
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Appendix 

A Symmetries of operator W 

We briefly discuss the requirements on the operator W in 
PI for the fulfillement of known strong interactions sym- 
metries. If U(g) is the unitary representation onto Hilbert 
space of an element g belonging to a symmetry group: 

U(g)WU(g)- 1 = W (A.l) 

Thus, as f] depends, by definition, only on Casimir opera- 
tors: 

U^PyUig)- 1 =P V 

\U(g)h v )(U(g)h v \ = IM(M (A.2) 

hv hv 

Since U (g) is a one-to-one correspondance, the require- 
ment is met if, for any g, every \U{g)hy) is a multi- 
hadronic state of the cluster, i.e. a \h' v ) or a linear com- 
bination of them. This is obvious if \hy) are eigenvectors 
of £/(<?), which is the case for the U(l) groups associated 
with abelian additive charges, and quite straightforward 
for isospin SU(2) since Py is the projector identity as far as 
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the isospin degrees of freedom are concerned; also charge 
conjugation symmetry is trivially satisfied. 

The situation is rather different for space-time symme- 
tries. In this case, the translation, rotation and reflection 
operators transform the projector Py in the projector onto 
the translated, rotated or reflected cluster respectively; 
only if this object is the same as the starting one, symme- 
try is fulfilled. Therefore, angular momentum and parity 
are conserved only if the cluster is spherical in shape, while 
energy and momentum are not conserved because of the 
finite volume. 



B Decomposition of the Poincare group 
projector 

The general transformation of the extended Poincare 
group g z may be factorized as: 



which is the Eq. 



C Proof of Equations (HU) and (|25|> 

We shall prove Eq. I|18|l in non-relativistic quantum me- 
chanics. It is assumed that |k) is a complete set of states in 
a region A with volume V , with eigenfunctions V'k(r) and 
that the transformation from a to r polarization states is 
unitary. Thus: 

ra\kr)=lt k{r)U " T S ! !H (C.l) 



if r i A 

where U aT is the element of a unitary matrix. Hence 
]T|<pa|kT)| 2 =£(pa|kT)(kT|pa) 



g z = T(x)ZA = T(x)ZL fi (OR (B- 1 ) = V / d 3 r / dV (pa|ra)(m|kr)(kr|rV)(r' ( T|p f 7 

t J A J A 

i~\ 4- t^t fr \ 4~ i Kit A- Ir-i ft Iaiii> t Tf-\s-i-i- f-\ -v ■ rvt / I I IP " } ' 



where T(x) is a translation by the four- vector x, Z = I, n is 
cither the identity or the space inversion and A = Lft(£)R 
is a general orthocronous Lorentz transformation written 
as the product of a boost of hyperbolic angle £ along the 
space-like axis n and a rotation R depending on three Eu- 
ler angles. Thus Eq @ becomes: 

' ^ ^ / ^x J dAD»(J(x)ZA)?U(T(x)ZA) 



A J A 



gip-(r'-r) 



E / d3r / dV^-^^ k (r)^(r')|C/ CTr | 2 (C.2) 



where we have used the normalization of the states 
(pIp') = <5 3 (P — P')- Since the tp^ are a complete set of 
eigenfunctions in A: 



2 ^ (27 
z=i.n v 



^V^(rM(r')=5 3 (r-r') 



(C.3) 



1 



\ " 



lim v 



2 ^ (2tt) 4 
xU(T{x))U(Z)U(A) 



d 4 x / dAe 1 ^7r 2 D ,y (A)" 



(B.2) 



where z = if Z = I and z = 1 if Z = 11. In the above 
equation, by dA we meant the invariant normalized mea- 
sure of the Lorentz group, which can be written as |2~2"] : 

dA = dL n (£)dR = sinh 2 £d£^-dR (B.3) 

4tt 

dR being the well known invariant measure of SU(2) 
group. 

If the initial state \i) has vanishing momentum, i.e. 
P = (M, 0) , then the Lorentz transformation A must not 
involve any non-trivial boost transformation with £ =/= 
for the matrix element D"(l\) 1 * not to vanish. Therefore 
A reduces to the rotation R and we can write: 



thus, taking into account that U is unitary, Eq. I|C.2|) turns 
to: 

J2 l(P-|kr)| 2 = 7^3 / d 3 r / dV <5 3 (r —0 = 7^ 
^ (2tt) 3 J a J a (2tt) 3 

(C.4) 

Likewise, the Eq. 1)28(1 can be proved by calculating: 
^(picri|kT)(kr|p 2 CT 2 ) = 



QED 



V / d 3 r / dV ( Pl ai|rai)(rai|kr)(kr|rV 2 ) 

7 1 J A J A 



{r a 2 \P2V2) 



-y 

2 ^ (2 
z=i,n v 



J d 4 x(2J+l) J dRe lP x TT z 



k,r 



= £ / d3r / d3r ' 



,ip 2 r -lpi-r 



A J A 



1 



(2tt) 3 



■Mr)^(r')U aiT U: 



xD J (R) x x *U(T(x))U(Z)U(R) (B.4) 

Since [Z, R] = 0, we can move the U{Z) operator to the 
right of U(R) and recast above equation as: 

1 



— x : / H 3 r P lr (P2-Pi) 

- °<ti><t 3 I u 1 



(C.5) 



P,J,X,1T 



(27T) 



d 4 x e^- x U{J{x)) 



:(2J+1) J dR D J {R) X X *U{R) l + ^ (n) ( B ,5) 



where the Eq. (|C.3(1 and unitarity of U have been used. 



D Extra strangeness suppression 

The use of an extra strangeness suppression parameter 75 
is quite common in statistical model analyses in canonical 
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and grand-canonical ensembles. We show here how this 
parameter can be inserted in the microcanonical ensem- 
ble giving rise to the usual formulae in the large volume 
limit. All that is needed is to multiply W by an opera- 
tor which add a factor 75 for each pair of valence strange 
quarks which is created or destroyed in the final state. 
Thus Eq. ® becomes: 



r 



f 



Nst\ 



(D.l) 



where N$i the number of strange quarks in the initial 
state and N$f — J2j ^j s i that in the final state, Sj being 
the number of valence strange quarks in the j th hadron 
species. Then, it is quite straightforward to extend the 
formulae shown in this paper for the presence of this ad- 
ditional factor. In particular, if Nsi = 0, the rate can be 
written: 



r 



where 



{Nj} 



(D.2) 
(D.3) 



and ^{Nj} as quoted throughout the paper. As far as 
mesons with fractional content Cs & [0, 1] of (ss) are con- 
cerned (77 for instance), an independent incoherent super- 
position of the rates is assumed as though the meson was a 
(ss) state in a fraction Cs of observed reactions. Therefore 
Eq. I|D.3(I must be rewritten as: 



(D.4) 



where 
fi = 



1 - 

s 3 

Is 



a 



Sj 



a 



Sjls 



for unflavoured mesons 
otherwise 

(D-5) 

These modifications lead to a corresponding modification 
of the microcanonical partition function: 



^' — / ] ^'{N 3 } - 
{Nj} 



{N 3 } 



(D.6) 



Under the same condition of validity of the approxima- 
tions, it can be proved, going along the equations quoted in 
Sect. 4, that this expression can be calculated explicitely. 
The grand-microcanonical partition function reads: 



fl = 



1 



x exp 



(27T) 4 

[ 



+ OC — IE 



dye 



(2 J, +1)V 



lP-y 



3 



(2tt) 3 



d 3 p log(l ± f je -^) ±l (D.7) 



The canonical partition function can be obtained start- 
ing from the above microcanonical partition function like 
in Sect. 5: 



Z(P, Q) 



1 



x exp 



3 



(2tt) m 
(2tt)3 



d 3 p log(l ± / je -/»-P-«u-*)±i 
(D.E 



which is the same as usually employed to derive the hadron 
multiplicities pp. 
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